Percolation with excluded small clusters and Coulomb blockade in a granular system 
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We consider dc-conductivity a of a mixture of small conducting and insulating grains slightly 
below the percolation threshold, where finite clusters of conducting grains are characterized by 
a wide spectrum of sizes. The charge transport is controlled by tunneling of carriers between 
neighboring conducting clusters via short "links" consisting of one insulating grain. Upon lowering 
temperature small clusters (up to some T-dependent size) become Coulomb blockaded, and are 
avoided, if possible, by relevant hopping paths. We introduce a relevant percolational problem of 
next-nearest-neighbors (NNN) conductivity with excluded small clusters and demonstrate (both 
numerically and analytically) that a decreases as power law of the size of excluded clusters. As a 
physical consequence, the conductivity is a power-law function of temperature in a wide intermediate 
temperature range. We express the corresponding index through known critical indices of the 
percolation theory and confirm this relation numerically. 

PACS numbers: 72.23.Hk, 73.22.-f, 72.80.Tm 
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A mixture of metallic and insulating grains is an im- 
portant practical system appearing in numerous appli- 
cations, such as ceramics, catalysts, and powders. It is 
also a useful model, describing many realistic composite 
materials (see [l| ) . Normally the properties of this sys- 
tem are described on the classical level: by the effective 
medium theory 0] (away from thepercolation threshold) , 
or by the percolation theory 0, H| (near the threshold). 
For traditional "old-fashioned" mixtures (with the sizes 
of grains a on the scale of microns), the classical approach 
is sufficient. However, in recent years a special attention 
was drawn to nanocomposite materials with character- 
istic a ~ 10 nm For such small grains essentially 
quantum effects should become important, especially on 
the insulating side of the percolation transition: tunnel- 
ing through insulating grains and the Coulomb blockade 
effect (see, e.g., @). 

In this paper, for definiteness, we consider a mixture of 
spherical grains of the same diameter a, the fraction x of 
them being metallic, while the fraction 1 — x being insu- 
lating (FigfTJ). Every pair of adjacent conducting grains 
is supposed to establish a good electric contact with di- 
mensionless conductance G 1, but still the resistance 
1/G of such a contact is larger than the resistance of a 
conducting grain itself, so that the latter will be alto- 
gether neglected in what follows. 

If one ignores the transport through insulating grains, 
then the macroscopic dc-conductivity a of the system 
identically vanishes below the percolation threshold (for 
x < x c ), so that the system undergoes a sharp metal- 
insulator transition (see 0, The infinite metallic 
cluster penetrating the entire system and carrying the 
dc current does not exist below the threshold: there are 
only finite metallic clusters embedded into the insulating 
matrix. The value of x c depends on the particular ar- 
rangement of the mixture. For example, for the regular 
densely packed 3d array of identical spheres (see Figfl] 




FIG. 1: A mixture of conducting grains (shown black) and 
insulating ones (white). Left panel: regular dense packing; 
right panel: topologically disordered random packing. 



left panel) x c w 0.2 fj], while for the 3d random pack- 
ing (right panel) x c ~ 0.35 0]. The critical exponents 
corresponding to the percolation transition are, however, 
universal and depend only on the dimensionality of the 
system d. 

Thus, to obtain a finite conductivity for x < x c one 
should necessarily take into account processes in which 
electrons traverse insulating grains. Let g be the dimen- 
sionless conductance of the shortest possible insulating 
bridge (consisting of a single insulating grain separating 
two conducting clusters, see Figl^). In what follows, 
we call these shortest bridges "links". Though usually 
g is exponentially small, its finiteness smears the metal- 
insulator transition, and the macroscopic conductivity of 
the system below the threshold becomes nonzero. 

If insulating grains are large and/or the temperature 
is relatively high, then electrons penetrate insulating 
grains due to activation mechanism (so that g = <? ac t oc 
exp{— W/T}, where W is the energy barrier for activa- 
tion of a carrier in the insulator). For this incoherent 
mechanism the conductance g^ of a bridge consisting 
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FIG. 2: Examples of paths (shown grey) through the insu- 
lator, connecting two large metallic clusters, (a) A single- 
grain bridge (a link). The links remain effective connectors 
of metallic clusters even at lowest temperatures; (b) A two- 
grain bridge. In the tunneling mode these and longer bridges 
do not contribute to the conductivity even at T > Ec; (c) 
A two-link chain path with an intermediate small (one-grain) 
metallic cluster. The paths involving small clusters are the 
first to be frozen out at T < Ec- 

of n insulating grains (see Fig^b for a n — 2 bridge) is 
g (n) = ffact/n, (1) 

and all such bridges with not very large n give compa- 
rable contributions to the effective conductance between 
two metallic clusters. Under the classical activation con- 
ditions the system is adequately described by a network 
made of randomly distributed classical conductances of 
two types: G 3> 1 and g <C 1, and the percolation 
through only large conductances is established at x = x c . 
The macroscopic conductivity er (measured in the units of 
(e 2 /h)a 2 ~ d ) of such a random network near the threshold 
is described by the formulas 0, Q 

~ 9{*c - x)-, A<°) « x c - x « 1, (2) 
cr£2 t ~ G(x - x c y, A<°) < x - x c < 1, (3) 

where fi is the universal conductivity exponent and s is 
the universal dielectric permeability exponent (see Table 
U). Here 

A< 0) = (. 9 /G) 1/( ^ +s) (4) 

is the width of the crossover region where the contribu- 
tion of the processes involving the insulating grains be- 
comes comparable with the conductivity of the infinite 
cluster: 

a^(* c -A(°>W^ c + A(°>). (5) 

Within the crossover region 

4rL ~ g»/^+»)G a '^ +s \ \x c -x\< A<°>. (6) 

The formula ^ does not involve g, since the current 
manages to totally avoid links in this range of concentra- 
tions, so that the potential drops are distributed over the 
metallic infinite cluster, while the insulating component 



of the mixture does not play any role in dc-conductivity. 
Hence the classical formula ((3]) is valid for description of 
the mixture independent on the nature of the processes 
in the insulator. That is not the case for the formula 
$Z$i that describes the regime where the potential drops 
occur only on the insulating links between the metallic 
clusters, while the latter are essentially equipotential. 
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0.875 
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TABLE I: Numerical values of some critical exponents 



For small grains and/or low temperatures, however, 
the tunneling mechanism dominates the transport via 
insulating grains, so that the conductance of the link is 
g = gtun oc exp{— 2ko}, where k is the tunneling decre- 
ment of the electronic wave function in the insulator [Io| . 
In the tunneling mode the effective conductance of an in- 
sulating bridge of n grains is 

= A n g? ua , (7) 

where A n is a model-dependent coefficient (typically 
A n <C 1 for n > 1). An electron in this case has to 
tunnel through a sequence of n insulating grains without 
any intermediate stop, and the corresponding tunneling 
amplitude exponentially decays with the length of the 
bridge. 

Under the tunneling conditions the model of a network 
with two types of conductances is generally inapplicable. 
Instead of one smeared percolation transition at x — x c , 
a cascade of smeared transitions at x — a;'™' arises, li™' 
being the threshold for percolation in the model where 
not only nearest neighbors but also next nearest and 
farther — up to m-th nearest neighbors — are assumed to 
be connected. In the range xi" 1 1 ' < x < Xc the con- 
ductivity a ~ g[™2 is dominated by bridges of length m. 
The crossover between the adjacent ranges at x = Xc 
is described by the formulas similar to (0) with g[™^ 
standing for g and g t ^™ standing for G. Note that this 
cascade of smeared phase transitions has the same ori- 
gin as the cascade-like concentration dependence of the 
temperature T g {x) of the spin-glass transition in those 
dilute magnetic alloys where the antiferromagnetic ex- 
change interaction of localized magnetic ions exponen- 
tially decreases with distance 

Since we are interested in the case of very small grains 
and low temperatures, in the present paper we will adopt 
the tunneling scenario and concentrate on the behavior 
of a just below the principal threshold 
the insulator-dominated regime where the direct nearest- 
neighbor (NN) percolation via metal grains is impos- 
sible but the next-nearest-neighbor (NNN) percolation 
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(which allows paths involving links, solitary insulating 
grains) still exists. Even without taking the Coulomb 
effects into account, the conductivity er ins in the tun- 
neling case should be somewhat suppressed compared 
to the classic percolational result (J2J): erj ns < cr^. In- 
deed, in the tunneling regime the contributions to the 
conductance between metallic clusters coming from the 
insulating bridges with n > 1 are suppressed, and only 
the contribution of links survives. Numerical simulations, 
however, demonstrate that a{ ns differs from cr^ only in a 
somewhat smaller numerical coefficient (see below) , while 
the index s remains the same: 

a inB (x) = ca>2 ~ g{x c ~ x)- s , c < 1. (8) 

This means that the problem with tunneling conduc- 
tances remains in the same universality class as the stan- 
dard one. The insensitivity of the index s to suppression 
of insulating long bridges indicates that the intercluster 
conductance in the case of activation scenario is dom- 
inated by relatively short bridges of n ~ 1 insulating 
grains anyway, while long bridges of length n 3> 1 are 
irrelevant. Note that only the contribution of these very 
long bridges could be responsible for the difference in the 
values of s for tunneling and activation cases (if there 
were any). 

Slightly below the percolation threshold the system 
consists of finite metallic clusters characterized by a wide 
spectrum of numbers of grains nt, sizes 7*j, conductances 
Qi 1 and charging energies Ei. The distribution function 
for clusters with n grains is well known (see, e.g., |4|): 



N(n) ~ n- T f(n/n c ), 



f(x) ~ e-«, (9) 



where df is the fractal dimension of the infinite cluster, 
and £ ~ a{x c — x)~ u 3> 1 is the correlation length. The 
clusters with the size r ~ £ and n ~ n c are commonly 
called critical clusters. The majority of grains belong to 
small clusters with n ~ 1 sites. The necessary exponents 
are given in the Table |TJ 

In the range 1 <C n <C n c ~ £ df the clusters are frac- 
tals, their physical characteristics are strongly fluctuat- 
ing. However, one can estimate the corresponding aver- 
age values for a fixed cluster size n. 

An average dimcnsionlcss conductance across the clus- 
ter i of fixed size n in the range 1 <C n <C n c can be 
estimated using standard expressions for conductivity of 
fractal objects (see Q). The result is: 

G(n) ee (£> n=const ~ Gn-Z/ d f , C = 2 - d + fi/v. 

The clusters with m < n ex t, where 

"cxt ~ G d f^ » 1, (10) 

typically have large conductances Qi > 1 and can be 
treated as structureless point-like objects. Clusters with 



rii > n ox t (and with Qi < 1) are, however, essentially 
extended objects. In this paper we consider only the case 



(11) 



when practically all existing clusters are point-like ones. 
The opposite case will be considered in a separate publi- 
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cation 

An average charging energy of clusters with fixed n 
can be estimated from the scaling law for the average 
screened Coulomb interaction 



(U(r)) 



E { °\a/ri v+s ^ v , 



where 



E, 



(0) 

c 



7" 



(12) 



(13) 



is the charging energy of a smallest one-grain cluster. 
The expression (|12j) is valid in the intermediate range 
1 <S r <C £ and gives a correct interpolation between 
the unscreened interaction U (r) ~ E@ jr at r ~ 1 and 
the screened one (U(r)) ~ E^/re at r > ^ (here I ~ 
(x c — x)~ s is the static long-range dielectric constant of 
the system [1, H3])- As a result, 



E c {n) ee (Ei) n=c 



(14) 



For temperatures T > Eq the Coulomb blockade ef- 
fects are completely suppressed, and the expression (JHJ, 
very similar to the classic result ([2]), is valid. 

In the range E^ < T < many clusters (those 
with Ei > T) are already frozen. The freezing condition 
may be rewritten as n < n(T), where 



n(T) ~ (4 0) /T) 



\vdf /(v+s) 



(15) 



How does the conduction proceed in such a system? 
Clearly, an electron would prefer to travel via large clus- 
ters and to avoid small frozen clusters. Is the NNN- 
percolation still possible after the exclusion of small clus- 
ters with n < n? How does the exclusion affect the con- 
ductivity (7? 

Let us call two metallic clusters "NNN-connected" if 
there is at list one link connecting these clusters (or, 
equivalently, at least one grain of the first cluster has 
a next nearest neighbor belonging to the second clus- 
ter). Now fix a certain n and consider the "refined" sys- 
tem of all finite clusters with "large" sizes n > n (see 
FigEbjC). If n ^> £ di is too large, then only a few of 
such clusters will be NNN-connected, and there will be 
no "NNN-percolation" through the system. On the other 
hand, for n ~ 1 the NNN-percolation should definitely 
exist. We define n c to be the value of h at which the 
NNN-percolation is first destroyed in the course of the 
refinement process. It is natural to expect 



(16) 
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FIG. 3: An infinite network of NNN-connected clusters. The 
smallest clusters (consisting of one grain) are shown light grey; 
the second smallest ones (of two grains) are shown dark grey; 
all larger metallic clusters are shown black, (a) The full sys- 
tem, (b) The system with excluded smallest clusters, (c) The 
system with excluded both smallest and next smallest clus- 
ters. Note that NNN-percolation exists in all three cases. 



On one hand, it is clear that n c can not be much larger 
than n c , since the distance between so extraordinarily- 
large clusters is exponentially large and they can not be 
NNN-connected. On the other hand, below we demon- 
strate that at n c ~ n c the NNN-percolation still exists. 
Indeed, the average intercluster distance R in the system 
of clusters with n in the range (1 — A)n c < n < n c : 



R~[N{n c )kn c } l l d ~kV d ii, 



(17) 



is of the same order of magnitude as the size of critical 
clusters £. (For derivation of (fT5|) we used the scaling 
relation r — 1 = d/dt between the critical exponents.) 
For A > A perc ~ 1 there is an infinite network of in- 
terpenetrating clusters (interpenetration of two clusters 
i and j can be defined, say, as a condition Rij < £). If 
A - A 



pcrc 



is not very small, then the fraction of critical 
clusters belonging to this interpenetrating network is of 
order unity. 

To estimate the average number Minks (£) of links 
NNN-connecting a typical pair of interpenetrating crit- 
ical clusters, we use the following reasoning (very close 
to the arguments of Coniglio [l3[ used for determina- 
tion of the fractal dimension of "red bonds"). Let us 
randomly choose some of the insulating sites (at concen- 
tration Ax > x c — x) and substitute conducting ones 
for them. Then the sample will be converted to a NN- 
percolating one with a similar value of correlation radius 
£. The initial critical clusters will become NN-connected 
with probability of order unity and form the infinite 
cluster in this new system. On the other hand, typi- 
cally two adjacent critical clusters become NN-connected 
when some insulating grain that was a critical link be- 
tween these clusters is replaced by a conducting grain. 
The probability of such an event is ~ Ax ■ Minks (£)■ Of 
course, one can also imagine a situation when the NN- 
connection between these two clusters is established via 
a "chain bridge" involving one or more smaller clusters, 
but the corresponding probabilities are proportional to 
higher powers of Ax and can be neglected. So, we con- 
clude that the NN-connection of adjacent critical clusters 
is typically established when Ax ■ Minks (£) ~ 1 and, con- 
sequently, 



M 



links 



(0 ~ l/(x c - x) » 1. 



(18) 



Thus, the interpenetrating critical clusters are typi- 
cally multiply NNN-connected, and therefore the infi- 
nite network of interpenetrating critical clusters is NNN- 
connected. 

Thus, if we exclude all clusters with n < n(T), the 
NNN-percolation in the system of remaining clusters per- 
sists if n(T) < n Cl or, equivalently, if T > Eq T \ 

Although the conductivity a of a "refined" system of 
clusters with sizes n > h is nonzero for n < n c (since the 
NNN-percolation is not destroyed), the magnitude of a is, 
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of course, suppressed by the refinement. For 1 <C h < n c 
this suppression is strong: 

a ins (x,n) ~al2(x)n- e . (19) 

The new critical index can be related to the known 
ones: 

(d-2)v + s-l 

e = ± '— . 20 

vdf 

To derive the relation (f2"0"|) we note that at h ~ n c our 
system consists of critical clusters only, each pair of adja- 
cent clusters being NNN-connected by ~ Minks (0 links, 
so that we can estimate 

cr ins (x, fi ~ n c ) ~ gN nnks {0£ 2 - d . (21) 

Then, comparing (]21|) with (| 19[) at n ~ n c , we imme- 
diately arrive at the relation (|20p . In two dimensions, 
making use of /i = s and Grassberger's precise result [3| 
/i/t/ = 0.9826(8), we obtain = 0.1227(4). Similarly, in 
3D one has « 0.28, which would be interesting to ver- 
ify numerically. The index increases with d and reaches 
= 1/2 at critical dimension d = 6. 

To verify the assumptions used in the above qualitative 
derivation, we performed moderate-statistics numerical 
computations for an ensemble of finite critical 2d sys- 
tems (at x = x c ) in the bus-bar geometry (see FigHJ. 
The sites adjacent to the two opposite sides of any sam- 
ple of the ensemble belong to the conducting subsystem 
by definition, while other sites are conducting with prob- 
ability x c and insulating with probability 1 — x c . All the 
conducting clusters with n < n are removed, but the bar- 
clusters, connected to the contacts, are "protected" — 
they are never removed. The current is allowed to flow 
from the left bar to the right one via large clusters with 
n > n, eventually jumping from cluster to cluster through 
the links. One half of these systems (those with the two 
bars directly connected to each other, so that the cur- 
rent actually avoids links) is percolating and has very 
large conductance oc G. The non-percolating systems of 
the other half have conductances oc g dominated by links. 
The expression for the size-dependent conductivity <7 cr , 
averaged over the non-percolating half of ensemble, is 

a cr (L,h) ~ gL s/v fC 6 ', Kn<L d f, (22) 

consistent with (|19p . The refinement process typically 
stops at n ~ L df , when the size of maximal removed 
cluster reaches the systems size L and only the two pro- 
tected bar-clusters (NNN-connected by ~ Mink(i) links) 
are left in the system (see FigHJj). 

Averages of finite conductivities <7 cr (L,n) for a L x L 
square vs. cutoff parameter n are given in FigO The 
data shown correspond to L = 256 for bond problem and 
to L = 512 for site problem. We used the Frank-Lobb 



(a) (b) 




FIG. 4: A finite critical non-percolating system in the bus- 
bar geometry, (a) Before any refinement: all NNN-connected 
clusters and corresponding links (white circles) are shown, 
(b) After ultimate refinement: only the two bar-clusters and 
critical links between them are left. 

algorithm to collect statistics over about 2 • 10 4 sam- 
ples per plot point. It is obvious that the deviations from 
power law Q22|) are small even for n ~ 1 or h adjacent to 
the saturation plateau; standard fits yield — 0.128(12) 
for the site variant and — 0.126(10) for the bond vari- 
ant, in good agreement with each other and with the 
theoretical prediction. 

We also checked numerically that NNN percolation 
indeed shares the value s with the ordinary NN prob- 
lem. Scaling of a c ^{L,n) with sample size L for L = 
16, 32 . . . 1024 suggested s/v = 0.983(9) in perfect agree- 
ment with percolation value s/v — 0.9826(8). 

The very fact that 0^0, reliably confirmed above both 
analytically and numerically, already gives us an impor- 
tant insight into the physics of the NNN-conduction in 
an infinite system with large but finite £. A priori one 
could imagine two possible scenarios of the conduction 
process: 

1 . The conducting paths in a non- refined sample con- 
sist predominantly of largest critical clusters with 
size n ~ n Cl directly connected by critical links. 
Smaller clusters with n -C n c are irrelevant. Were 
this scenario the actual one, the refinement pro- 
cess would not lead to any considerable suppres- 
sion of the conductivity as long as n C n c , and, 
consequently, one would expect 9 = and s = 
1 - (d-2)u. 

2. Though the number of relevant conducting paths in 
a non-refined sample depends on x and diverges as 
x — > x c , each path consists predominantly of small 
clusters with size n ~ 1. When all clusters with n < 
h are removed from our system due to refinement, 
the leading role passes on to the smallest remaining 
clusters with n ~ n. For this scenario one expects 
9 > and s > 1 - (d-2)v. 

Our results definitely indicate that the second scenario, 
and not the first one, is real. Note that the opposition 
of the two concepts described above is reminiscent of the 
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L=256, bond x 
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FIG. 5: Average conductivity of a L x L square as a function 
of the cutoff parameter fi for site and bond models. Lines 
correspond to fits with 6 — 0.128 (bond) and 6 — 0.126 (site). 



alternative between the oversimplified "links and nodes" 
model 16] and a more sophisticated "links and blobs" 
model [17[ in the standard NN-percolation problem . 

Thus, substituting (fl~5j) into (fT9|) . for the temperature 
dependence of the conductivity we finally obtain 

a ins (x,T) ~ a ins (z,n(T)) ~ a^l(x)(T/E^f , (23) 



e = 



d f v9 _ (d-2)u + s-l 



v + s v + s 

valid in the wide intermediate temperature range 



4 cr) «t«4 0) , 



(24) 



(25) 



and not very close to the percolation transition (outside 
the critical crossover domain): 

A(T) < x c -x < 1, A(T) = A<°) • (T/E { ° ] ) e /^+ s l 



Within the critical crossover region (at \x c 
the conductivity stops to depend on x: 

& = Qn/(p + s). 



< 



A(T)) 



(26) 



the 



The phase diagram for different regimes in 
temperature-concentration plane is shown in FigJS] 

At temperatures T < E y c the T-dependence of the 
conductivity becomes stretched exponential; it is due to 
the fact that at these low temperatures n(T) > n c and 
the current has no chance to avoid the Coulomb block- 
aded clusters. The underlying physics is very different 
from the exclusion of the clusters considered above; it is 
discussed in a separate publication 
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FIG. 6: Domains of the power-law T-dependence on the 8 — t 
plane. Here 8 = x c — x is the distance from the percola- 
tion transition, t = T/Eq is the dimensionless temperature, 
text ~ G ("+»)/[(2-d)^+M] < L The upper crossover li ne is 
given by 8 ~ t 2 /("+ 3 ), the lower one— by 5 ~ A (0) f e/(fl+s) . 



In conclusion, we have demonstrated that the Coulomb 
blockade effect in a granular material close to the per- 
colation threshold leads to the power-law temperature 
dependence of conductivity, valid in a wide intermedi- 
ate temperature-range. This dependence arises due to 
the wide power-law distribution of the charging energies 
of conducting clusters, characteristic for the critical re- 
gion near the percolation threshold. The corresponding 
exponent is universal, it depends only on the dimension- 
ality of the system and can be expressed through known 
critical exponents of the percolation theory. It is inter- 
esting to note that in three dimensions both exponents 
6 = 0.38 and 6' = 0.28 are close to the exponent 1/3 
used for the fitting of the experimentally observed critical 
T-dependence of conductivity of thick amorphous (but, 
presumably, not granulated) films of InO and TiN, in 
which the superconductivity was suppressed by strong 
magnetic field [3]. It is possible that for the latter sys- 
tems the similar power law cr(T)-dependence is due to the 
self-induced inhomogeneity of disordered systems devel- 
oping at low temperatures when the system approaches 
the MIT (see 0). 

The Authors are indebted to M.V.Feigelman and 
V.F.Gantmakher for valuable discussions of physics of 
granular metals. This work was supported by RFBR 
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